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Axion Like Particles (ALPs), postulated to solve the strong-CP problem, are predicted to couple 
with photons in the presence of magnetic fields, which may lead to a significant change in the 
observed spectra of gamma-ray sources such as AGNs. Here we simultaneously consider in the same 
framework both the photon/axion mixing that takes place in the gamma-ray source and that one 
expected to occur in the intergalactic magnetic fields. We show that photon/axion mixing could 
explain recent puzzles regarding the observed spectra of distant gamma-ray sources as well as the 
recently published lower limit to the EBL intensity. We finally summarize the different signatures 
expected and discuss the best strategy to search for ALPs with the Fermi satellite and current 
Cherenkov telescopes like CANGAROO, HESS, MAGIC and VERITAS. 



I. PHOTON/AXION OSCILLATIONS 



Axions were postulated to solve the strong-CP 
problem in QCD in the 1970s 0, and are still the 
most compelling explanation for it. Moreover, they 
are also appealing for astrophysicists, since they are 
valid dark matter candidates to constitute a portion or 
the totality of the non-barionic cold dark matter con- 
tent predicted to exist in the Universe. But probably 
the most interesting property of axions, or in a more 
generic way, Axion-Like Particles (ALPs), where the 
mass TOa and the coupling constant are not related to 
each other, is that they are expected to oscillate into 
photons (and viceversa) in the presence of magnetic 
fields [1, 0] ■ This oscillation of photons to ALPs are 
the main vehicle used at present in axion searches, like 
those carried out by CAST ^ or ADMX ^, but they 
could also have important implications for astronom- 
ical observations. For example, they could distort the 
spectra of gamma-ray sources, such as Active Galactic 
Nuclei JAlGNs) 0-[i3 or galactic sources in the TeV 
range [ll| . These distortions may be detected by cur- 
rent gamma-ray experiments, such as Imaging Atmo- 
spheric Cherenkov Telescopes (lACTs) hke MAGIC 
ll, HESS [13 , VERITAS ;i4j or CANGAROO-III 
15| . covering energies in the range 0.1-20 TeV, and 
by the Fermi satellite [l3|, which operates at energies 



in the range 0.02-300 GeV. 

Given a domain of length s, where there is a roughly 
constant magnetic field and plasma frequency, the 
probability of a photon of energy Ej to be converted 
into an ALP after traveling through it can be written 
as 



2 sin'(Ao.c s/2) 

Here Aosc is the oscillation wave number: 
Al^^ ~ {AcM + Api - Aaf + 4A 



(1) 



(2) 



Ab that gives us an idea of how effective is the mixing, 
i.e. 



2 M 



1.7 X 10-2^ Mil B„,G cm-^, (3) 



where Bt the magnetic field component along the po- 
larization vector of the photon and Mn the inverse of 
the coupling constant. 

AcM is the vacuum Cotton-Mouton term, i.e. 



-1.3 X 10-^' Bla cnir\ (4) 



* A summary of the work presented in Ref . Q 



where Bcr = mg/e ~ 4.41 x 10^'^ G the critical mag- 
netic field strength (e is the electron charge). 
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Api is the plasma term: 

(5) 

where Wpi = yj Airane/rne = 0.37 x 

jieV \J ne/cm~^ the plasma frequency, 
the electron mass and the electron density. 
Finally, Aa is the ALP mass term: 

A, = ^ ^ 2.5 X lO-^Vl,,^ cm-\ (6) 

Note that in Eqs.(I5][ni) we have introduced the di- 
mensionless quantities BmG = B/10~^ G, AIu = 
M/10" GeV and m^ev = m/lQ-^ eV. 

Since we expect to have not only one coherence do- 
main but several domains with magnetic fields differ- 
ent from zero and subsequently with a potential pho- 
ton/axion mixing in each of them, we can derive a 
total conversion probability [TT'| as follows: 

P^^a^l[l-eM-iNPo/2)] (7) 

where Pq is given by Eq.([T|) and N represents the 
number of domains. Note that in the limit where 
N Pq —i' oo, the total probability saturates to 1/3, 
i.e. one third of the photons will convert into ALPs. 

It is useful here to rewrite Eq. ((T|) following Ref . Q , 
i.e. 



where Spc = s/parsec. Some astrophysical environ- 
ments fulfill the above mixing requirement and the 
Mil constraints imposed by CAST. Indeed, when us- 
ing A/ii > 0.114 in Eq. ([TT|), we can deduce that as- 
trophysical sources with • Spc > 0.01 will be valid. 
This product B-s also determines the maximum en- 
ergy Eimax to which sources can accelerate cosmic rays, 
and is given by Emax = 9.3 x 10^° Bq ■ Spc eV (Hillas 
criterium). Since we observe cosmic rays up to 3 x 
10^° eV, Bg Spc can be as high as 0.3, which also 
means that sources with Bg-Spc =0.01 are completely 
plausible and should exist. Therefore, photon/axion 
mixing may have important implications for astro- 
nomical observations (AGNs, pulsars, GRBs...). Note, 
however, that an efficient mixing can be expected to 
occur not only in compact sources: the mixing will 
be also present in the Intergalactic Magnetic Fields 
(IGMFs), with typical values of '-^l nG for the B field, 
provided that the source is located at cosmological dis- 
tances (Spc = 10*) in order to ensure that B -s is still 
larger than 0.01. 

Therefore, in order to correctly evaluate the pho- 
ton/axion mixing effect for distant sources, it will be 
necessary to handle under the same consistent frame- 
work the mixing that takes place inside or near the 
gamma-ray sources together with that one expected 
to occur in the IGMF. In the literature, both effects 
have been considered separately. We neglect, how- 
ever, the mixing that may happen inside the Milky 
Way due to galactic magnetic fields. [29| 



Po = 



sm 




(8) 

so that we can define a characteristic energy, Ecrit, 
given by: 



Er. 



M 



2 B 



or in more convenient units: 



Ecr^t{GeV) 



Mil 



0.4 Bg 



(9) 



(10) 



where the subindices refer again to dimensionless 
quantities: nif^eV = m/^eV, Mn = M/10^^ GeV 
and Bg = B/ Gauss; m is the effective ALP mass 
= \m1—ujpi\. Recent results from the CAST exper- 
iment [5| give a value of Mu > 0.114 for axion mass 
TUa < 0.02 eV. Although there are other limits derived 
with other methods or experiments, the CAST bound 
is the most general and stringent limit in the range 
10-" eV < TOa < 10-^ eV. 
In order to have an efficient conversion, we need 



15 Bg s. 



pc 



> 1 



(11) 



A. Mixing in the source 

To illustrate how the photon/axion mixing inside 
the source works, we show in Figure [T] an example for 
an AGN modeled by the parameters listed in Table HI 
The only difference is the use of an ALP mass of 1 
fj,eV instead of the value that appears in that Table, 
so that we can obtain a critical energy that lies in the 
GeV energy range; indeed, we get Ecrit =0.19 GeV 
according to Eq. (ITU|) . Note that the main effect just 
above Ecru is an attenuation of the total expected flux 
coming from the source. [30| 



B. Mixing in the IGMFs 

As already discussed, despite the low magnetic field 
B, the photon/axion oscillation can take place in the 
IGMFs due to the large distances. However, the 
strength of the IGMFs is expected to be many or- 
ders of magnitude weaker (^nC) than that of the 
source and its surroundings (~G). Consequently, as 
described by Eq. ([TU)) . the energy at which pho- 
ton/axion oscillation occurs in this case is many orders 
of magnitude larger than that at which oscillation can 
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FIG. 1: Example of photon/axion oscillations inside the 
source or vicinity, and its effect on the source intensity 
(solid line), which was normalized to 1 in the Figure. We 
used the parameters given in Table |I] to model the AGN 
source, but we adopted an ALP mass of 1 fieV. This gives 
Ecrit ~ 0.19 GeV. The dot-dashed line represents the max- 
imum (theoretical) attenuation, i.e. 1/3. 

occur in the source and its vicinity. Assuming a mid- 
value of B~0.1 nG, and Mn = 0.114 (CAST upper 
limit), the effect could be observationally detectable 
by current lACTs {Ecru < 1 TeV) only if the ALP 
mass is ma < 6 x 10""'^'' eV, i.e. we need ultra-light 
ALPs. For example, we get Ecrit = 28.5 GeV when 
using mjj = 10^^" eV, which is the value given in Ta- 
ble |T] as our reference one. 

TABLE L Parameters used to calculate the total pho- 
ton/axion conversion in both the source (for the two AGNs 
considered, 3C 279 and PKS 2155-304) and in the IGM. 
The values related to 3C 279 were obtained from Ref. [l3| , 
while those ones for PKS 2155-304 were obtained from 
Ref. [isj . This Table represents our fiducial model. 





Parameter 


3C 279 


PKS 2155-304 




B (G) 


1.5 


0.1 


Source 


ed (cm"^) 


25 


160 


parameters 


L domains (pc) 


0.003 


3 X 10"'' 




B region (pc) 


0.03 


0.003 




z 


0.536 


0.117 


Intergalactic 


ed,int (cm~^) 


10-^ 


10-^ 


parameters 


B^r,t (nG) 


0.1 


0.1 




L domains (Mpc) 


1 


1 


ALP 


M (GeV) 


1.14 X 10^" 


1.14 X 10^" 


parameters 


ALP mass (eV) 


io-i« 


10-10 



In our model, we assume that the photon beam 
propagates over N domains of a given length, the mod- 



ulus of the magnetic field roughly constant in each of 
them. We take, however, randomly chosen orienta- 
tions, which in practice is also equivalent to a variation 
in the strength of the component of the magnetic field 
involved in the photon/axion mixing Moreover, 
it will be necessary to include the effect of the Ex- 
tragalactic Background Light (EBL) in the equations 
as well, its main effect being an additional attenua- 
tion of the photon flux, especially at energies above 
100 GeV. Indeed, the EBL plays a crucial role to cor- 
rectly evaluate and understand the importance of the 
intergalactic mixing. The induced effect can be an 
attenuation or an enhancement of the photon flux at 
Earth, depending on distance, magnetic field and the 
EBL model considered (see next Section). 



C. Source and intergalactic mixings working 
together 

In conclusion, AGNs located at cosmological dis- 
tances will be affected by both mixing in the source 
and in the IGMF. As already mentioned, up to now 
previous works have focused only in studying the pho- 
ton/axion mixing either inside the source or in the 
IGMFs. Instead, for the first time we carried out a 
detailed study of the mixing in both regimes under the 
same consistent framework. In order to observe both 
effects in the gamma-ray band, we need ultra-light 
ALPs. That is the reason why in our fiducial model, 
presented in TablelH we use a value of m^, = 10"^° eV, 
which implies Ecrit ~ 30 GeV for the IGMF mixing 
(for B'^0.1 nG) and Ecrit ^ 1 eV within the source 
and its vicinity (B^l G). Consequently, both effects 
need to be taken into account. We show in Fig. [5] 
a diagram that outlines our formalism. Very sque- 
matically, the diagram shows the travel of a photon 
from the source to the Earth in a scenario with ALPs. 
In the same Figure, we show the main physical cases 
that one could identify inside our formalism: mixing 
in both the source and the IGMF, mixing in only one 
of these environments, the effect of the EBL, etc. 



II. AXION BOOSTS 

Since we expect the intergalactic mixing to be more 
important for larger distances, due to the more promi- 
nent role of the EBL, we chose two distant astrophysi- 
cal sources as our benchmark AGNs: the radio quasar 
3C 279 (z=0.536) and the BL Lac PKS 2155-304 at 
z=0.117. We summarize in Table|T]the parameters we 
have considered in order to calculate the total pho- 
ton/axion conversion in both the source and in the 
intergalactic medium. As for the EBL model, we chose 
the Primack [l9| and Kneiske best-fit [20| ones. They 
represent respectively one of the most transparent and 
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FIG. 2: Sketch of the formahsm here presented, where both mixing inside the source and mixing in the IGMF are 
considered under the same consistent framework. Photon to axion oscillations (or vice-versa) are represented by a 
crooked line, while the symbols 7 and a mean gamma-ray photons and axions respectively. This diagram collects the 
main physical scenarios that we might identify inside our formalism. Each of them are squematically represented by a 
line that goes from the source to the Earth. 



one of the most opaque models for gamma-rays, but 
still within the limits imposed by the observations. 




0.1 1.0 10.0 100.0 1000.0 

E (QeV) 



Primack B = 0,1 


nG PKS 2155-304 ■ 


Primock B-1 nG 




Primack B^O.Ol 


nG 
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t B=l nG / 
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B=0.01 nG .■■ /- 


1 1 -s^-- 
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FIG. 3: Boost in intensity due to ALPs for the Kneiske 
best-fit (red light lines) and Primack (blue dark lines) EBL 
models, and for different values of the magnetic field. We 
used those parameters presented in Table |T] for 3C 279 
(z=0.536) and PKS 2155-304 (z=0.117). 

In order to quantitatively study the effect of ALPs 



on the total photon intensity, we plot in Figure [3] the 
difference between the predicted arriving photon in- 
tensity without including ALPs and that one obtained 
when including the photon/axion oscillations (called 
here the axion boost factor). This was done for our 
fiducial model (Table IJ) and for the two EBL mod- 
els considered. The plots show differences in the ax- 
ion boost factors obtained for 3C 279 and PKS 2155- 
304 due mostly to the redshift difference. The in- 
ferred critical energies for the mixing in the source 
are Ecrit = 4.6 eV for 3C 279 and Ecr^t = 69 eV for 
PKS 2155-304, while for the mixing in the intergalac- 
tic medium we obtain Ecrit = 28.5 GeV (the same 
for both objects). For B=0.1 nG, and in the case of 
3C 279, the axion boost is an attenuation of about 
16% below the critical energy (due to mixing inside 
the source). Above this critical energy and below 200- 
300 GeV, where the EBL attenuation is still small, 
there is an extra attenuation of about 30% (mixing 
in the IGMF). Above 200-300 GeV the axion boost 
reaches very high values: at 1 TeV, a factor of ~7 
for the Primack EBL model and ^^340 for the Kneiske 
best-fit model. We find that the more attenuating 
the EBL model considered, the more relevant the ef- 
fect of photon/axion oscillations in the IGMF, since 
any ALP to photon reconversion might substantially 
enhance the intensity arriving at Earth. In the case 
of PKS 2155-304, the situation is different from that 
of 3C 279 due to the very low photon-attenuation at 
the source and, mostly, due to the smaller source dis- 
tance. Furthermore, a very interesting result is found 
when varying the modulus of the intergalactic mag- 
netic field. Higher B values do not necessarily trans- 
late into higher photon flux enhancements. There is 
always a B value that maximizes the axion boost fac- 
tors; this value is sensitive to the source distance, the 
considered energy and the EBL adopted model (see 
Ref. [ll for a detailed discussion on this issue). 
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lACTs observations 

Look for systematic intensity 
enhancements at energies where 
the EBL is important. 

Distant (z > 0.2) sources at the 

highest possible energies (>1 
TeV), to push EBL models to the 
extreme. 

Source and EBL model 
dependent, but very important 
enhancement expected in some 
cases. 



Fermi/IAT 

Look for intensity drops in the residuals 
("best-model"-data). 
Source model dependent. 
Powerful, relatively near AGNs. 



Fermi/UT and/or lACTs 

Look for intensity drops in the residuals. 

Only depends on the IGMF and axion properties (mass 

and coupling constant). 

independent of the sources -> CLEAR signature! 



FIG. 4: Summary of the best observational strategy to look for ALPs with Fermi and lACTs. See Ref. [2| for a deeper 
discussion. 



III. DETECTION PROSPECTS FOR FERMI 
AND lACTS 

If we accurately knew the intrinsic spectrum of the 
sources and/or the density of the EBL, we should be 
able to observationally detect axion signatures or to 
exclude some portions of the parameter space. We 
lack this knowledge, so detection is challenging but 
we believe that still possible. The combination of the 
Fermi/LAT instrument and the lACTs, which cover 
6 decades in energy (from 20 MeV to 20 TeV) is very 
well suited to study the photon/axion mixing effect. 
Nevertheless, and before assuming an scenario with 
axions to interpret the observational data, we should 
try to describe the observational data (preferably sev- 
eral AGNs located at different redshifts, as well as the 
same AGN undergoing different flaring states, from 
radio to TeV) with "conventional" theoretical models 
for the AGN emission and for the EBL. If these "con- 
ventional" models for the source emission and EBL 
fail (i.e. if we have important residuals for the best-fit 
model), then the axion scenario should be explored. 
Fig. m summarizes what could be a good observational 
strategy to look for ALPs with Fermi and lACTs. 



IV. ARE WE DETECTING AXIONS 
ALREADY? 

Recent gamma observations might already pose sub- 
stantial challenges to the conventional models to ex- 
plain the observed source spectra and/or EBL density. 



e.g.: 

• The VERITAS Collaboration recently claimed 
a detection above 0.1 TeV coming from 3C66A 
(z=0.444). The EBL-corrected spectrum seems 
to be harder than 1.5 [2l|. 

• TeV photons coming from 3C 66A? (see Refs. 
[23 - [24| ). If so, difficult to explain with conven- 
tional EBL models and physics. 

• The lower limit on the EBL at 3.6 /im was re- 
cently revised upwards by a factor ^2, suggest- 
ing a more opaque universe [25j . 

• Some sources at z = 0.1-0.2 seem to have harder 
intrinsic energy spectra than previously antici- 
pated Hi]. 

While it is still possible to explain these observa- 
tions with conventional physics, the axion/photon os- 
cillation would naturally explain these puzzles, since 
we get more high energy photons than expected as 
well as a softer intrinsic spectrum in a scenario with 
ALPs. An example is given in Fig. [5] (from Ref. (27j). 
where the effect of the existence of ALPs in the intrin- 
sic spectrum of 3C 279 is shown assuming the fiducial 
model presented in Table U and a Kneiske best-fit EBL 
model. 



V. CONCLUSIONS 

If axions exist, then they could distort the spec- 
tra of astrophysical sources importantly. In addition, 
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3C279, Kneiske+04 bestfit, B=0.1nG 
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FIG. 5: Effect of the existence of ALPs in the intrinsic 
spectrum of 3C 279, assuming the fiducial model given in 
Table [I] and a Kneiske best-fit model for the EBL. Blue 
points correspond to observational data taken with the 
MAGIC telescope [H]. 



if photon/axion mixing occurs in the IGMFs, then a 
mixing in the source will be at work as well. In par- 
ticular, for axion masses of the order or 10"^" eV, the 
induced effect is expected to be present in the gamma 
ray energy range. 

Since photon/axion mixing in both the source and 
the IGM are expected to be at work over several 
decades in energy, it is clear that a joint effort of Fermi 
and current lACTs is needed. The Fermi/LAT instru- 
ment is expected to play a key role, since it will detect 
thousands of AGNs (up to z~5), at energies where the 
EBL is not important. By the other hand, lACTs will 
be specially important at higher energies (>300 GeV), 
where the EBL is present. 

Main caveats: the effect of photon/axion oscillations 
could be attributed to conventional physics in the 
source and/or propagation of the gamma-rays towards 
the Earth. However, we believe that detailed observa- 
tions of AGNs at different redshifts and different flar- 
ing states could be successfully used to identify the 
signature of an effective photon/axion mixing. 



[1] Sanchez-Conde M. A., Paneque D., Bloom E., Prada 
F. and Domfnguez A., 2009, Phys. Rev. D, 79, 123511 

[2] Peccei R. D. and Quinn H. R., 1977, Phys. Rev. Lett, 
38, 1440 

[3] Dicus D. A., Kolb E. W., Teplitz V. L. and Wagoner 

R. v., 1978, Phys. Rev. D, 18, 1829 
[4] Sikivie P., 1983, Phys. Rev. Lett, 51, 1415 [Erratum 

ibid., 1984, Phys. Rev. Lett, 52, 695 
[5] Andriamonje S. et al, 2007, JCAP, 0704, 010 
[6] Duffy L. D. et al., 2006, Phys. Rev. D, 74, 012006 
[7] Hooper D. and Serpico P., 2007, Phys. Rev. Lett, 99, 

231102 

[8] De Angelis A., Roncadelh M. and Mansutti O., 2007, 

Phys. Rev. D, 76, 121301 
[9] Hochmuth K. A., Sigl G., 2007, Phys. Rev. D, 76, 

123011 

[10] Simet M., Hooper D. and Serpico P., 2008, Phys. Rev. 
D, 77, 063001 

[11] Mirizzi A., Raffelt G. G. and Serpico P., 2007, Phys. 

Rev. D, 76, 023001 
[12] Lorentz E. et al., 2004, New Astron. Rev., 48, 339 
[13] Hinton J. A., 2004, New Astron. Rev., 48, 331 
[14] Weekes T. C. et al, 2002, Astropart Phys., 17, 221 
[15] Enomoto R. et al., 2002, Astropart. Phys., 16, 235 
[16] Gehrels N. and Michelson P., 1999, Astropart Phys., 

11, 277 

[17] Hartmann et al, 2001, ApJ, 553, 683 
[18] Kusunose M. and Takahara P., 2008, ApJ, 682 
[19] Primack J.R., 2005, Proceedings of the Gamma 
2004 Symposium on High Energy Gamma Ray As- 
tronomy, 26-30 July 2004, Heidelberg, Germany, 



astro- ph/0502177| 
[20] Kneiske T. M., Bretz T., Mannheim K. and Hartmann 

D. H., 2004, A&A, 413, 807 
[21] Acciari V. A., 2009, Ap,J, 693, L104 
[22] Aliu E. et al., 2009, ApJ, 692, L29 
[23] Neshpor Y. I., Stepanyan A. A., Kalekin O. P., Fomin 

V. P., Chalenko N. N. & Shitov V. G. 1998, Astron. 

Lett, 24, 134 

[24] Stepanyan A. A., Neshpor Y. I., Andreeva N. A., 
Kalekin O. P., Zhogolev N. A., Fomin V. P.& Shitov 
V. G., 2002, Astron. Rep., 46, 634 

[25] Levenson L. R. & Wright E. L., 2008, Ap,J, 683, 585 

[26] Krennrich, Dwek and Imran, 2008, ApJ, 689, L93 

[27] Sanchez-Conde M. A. & Domfnguez A., in prep. 

[28] Albert et al., 2008, Science, 320, 1752 

[29] At present, a concise modeling of this effect is still 
very dependent on the largely unknown morphology 
of the magnetic field in the galaxy. Furthermore, in 
the most optimistic case, it would produce a photon 
flux enhancement of ~3% of the initial photon fiux 
emitted by the source [Tol ]. 

[30] Note, however, that the attenuation starts to decrease 
at higher energies (>10 GeV) gradually, the reason 
being the crucial role of the Cotton-Mouton term, 
which makes the efficiency of the source mixing to 
decrease as the energy increases. See Ref. [H for de- 
tails. 

[31] We refer to Ref. 1] for a detailed description of the 
model and the related equations. 



eConf C091122 



